Quantum trajectory phase transitions in the micromaser 
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We study the dynamics of the single atom maser, or micromaser, by means of the recently in- 
troduced method of thermodynamics of quantum jump trajectories. We find that the dynamics of 
the micromaser displays multiple space-time phase transitions, i.e., phase transitions in ensembles of 
quantum jump trajectories. This rich dynamical phase structure becomes apparent when trajecto- 
ries are classified by dynamical observables that quantify dynamical activity, such as the number of 
atoms that have changed state while traversing the cavity. The space-time transitions can be either 
first-order or continuous, and are controlled not just by standard parameters of the micromaser but 
also by non-equilibrium "counting" fields. We discuss how the dynamical phase behavior relates to 
the better known stationary state properties of the micromaser. 
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INTRODUCTION 



Dynamical transitions are a common feature of driven 
open quantum systems. Probably the best known ex- 
ample is the laser whose behavior close to threshold has 
many features in common with a thermodynamic sys- 
tem close to a continuous phase transition [T]. Another 
text-book example which has an especially rich dynami- 
cal phase behavior is the micromaser [IHS] where a flux 
of atoms pass one at a time through a microwave cav- 
ity. As the flux of atoms is increased the cavity un- 
dergoes a continuous dynamical crossover, followed by 
a series of flrst-order like crossovers. These crossovers 
are in practice rounded transitions that only become for- 
mally analogous to phase transitions in the limit where 
the number of atoms sent through the cavity during its 
lifetime becomes inflnitely large [3H5]. The dynamics of 
the micromaser has been investigated very successfully in 
experiment by measuring the state of the atoms emerg- 
ing from the cavity [3], an approach which gives access 
to quantum trajectories of the system [51[7]. 

In this work we investigate the connection between the 
dynamical transitions in the micromaser and thermody- 
namic phase transitions in a new way by analysing the 
statistical properties of quantum trajectories of the cav- 
ity held. We build on a recent Letter [S] where two of 
us introduced a thermodynamic formalism for the study 
of quantum jump trajectories [7] of open quantum sys- 
tems [5] . This formalism |S] , based on the large-deviation 
(LD) method |TD], treats statistical ensembles of dynam- 
ical trajectories in a manner equivalent to the way in 
which equilibrium statistical mechanics treats ensembles 
of conflgurations. This method allows the statistics of 
dynamical events to be understood using the conceptual 
framework of thermodynamics [11' . In the classical con- 
text, this "statistical mechanics of trajectories" approach 
has proved useful in uncovering space-time phase transi- 
tions (i.e. transitions in the space of temporal trajecto- 
ries) in the dynamics of slowly relaxing systems such as 
glasses [12] ■ In Ref. [8] we discussed how similar space- 
time transitions can be present in dissipative quantum 
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1: (Color online) Space-time phase diagram of the mi- 
cromaser. The dynamics displays multiple space-time transi- 
tions between dynamical phases of distinct dynamical activity 
k. The phase diagram is in terms of the "pump parameter" 
a, the standard control parameter of the micromaser, and 
the "counting field" s (see text for details). The white lines 
indicate first-order dynamic phase boundaries. The right- 
most phase boundary ends at a critical point with Sc > 
and ac> 1. 



systems, and showed as an example one such transition 
in the micromaser. Here we go beyond this result by 
providing a comprehensive study of the dynamical phase 
structure of the micromaser using the LD approach. 

Figure [l] summarizes the findings we describe in detail 
below. It shows the "space-time phase diagram" of the 
micromaser obtained via the method of [BJ . The figure is 
a density map for the average number of quantum jumps 
[7] corresponding to the average number of atoms, (fc), 
which have released a quantum of energy into the cavity, 
per unit time and at stationary conditions. This quan- 
tity is a dynamical order parameter, as it allows distinct 
dynamical phases to be classified and distinguished be- 
tween. The figure shows the dependence of (k) on two pa- 
rameters, the reduced atom pump rate a (defined below) 
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and the non-equilibrium "counting" field s f5|. Whilst 
a can be varied directly in an experiment, the behav- 
ior at different values of the parameter s describes the 
statistics of quantum jumps whose probability is biased 
by the value of s so that they become more (s < 0) or 
less (s > 0) likely to occur. The behavior at s = de- 
scribes the statistics of the trajectories which would be 
observed in an experiment. Remarkably, we find that at 
certain values of a and s the order parameter changes in 
a way which appears to be singular^ even though in this 
case only 100 atoms go through the cavity in its lifetime. 
These points, which occur at non-zero values of s, are the 
location of phase transitions between phases of distinct 
"dynamical activity" [8l[l2Hl4]. 

The rest of this paper is organized as follows. In Sec. 
II we review the standard theoretical description of the 
micromaser and describe the quantum trajectories gener- 
ated by measuring the state of the emerging atoms. We 
then introduce the counting field and describe how large 
deviation theory can be used to calculate the statistical 
properties of the trajectories. We use a mean-field ap- 
proach to uncover the statistical properties of the trajec- 
tories for arbitrary s in Sec. Ill, and hence obtain phase 
diagrams like that shown in Fig. [T] Then in Sec. IV we 
compare the mean-field results with an exact numerical 
analysis. Finally, Sec. V contains our conclusions and a 
short discussion. 
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FIG. 2: (a) Schematic operation of the micromaser The 
cavity mode interacts with a thermal bath, and is driven by 
two level atoms pumped, one by one, in their excited state. 
When the atoms traverse the cavity they can decay to their 
ground state, releasing a photon that is absorbed by the cav- 
ity, (b) Mean photon number (n) in the cavity in the steady 
state (in units of A'ex) as a function of the "pump parameter", 
a = (jiy^N^. Here iVcx = 100 and u = 0.15. (c) The Fano 
factor. 



II. MODEL AND FORMALISM 

The micromaser consists of a microwave cavity which is 
pumped by a sequence of two level atoms which are sent 
into the cavity at a constant (pump) rate; the cavity is 
also coupled to a thermal bath [U [SHS] . In the simplest 
case, which we will consider here, the atoms are all pre- 
pared in the excited state and pass through the cavity 
one at a time. This setup is sketched in Fig. |2ja). The 
atoms are assumed to interact resonantly with a single 
mode of the cavity. The cavity mode reaches a steady 
state which is sensitive to the pump rate and the atom- 
cavity coupling as well as the properties of the thermal 
bath. 

A standard treatment [H [5] for the micromaser focuses 
on the dynamics of the cavity. After integrating out the 
thermal bath and the atom degrees of freedom, and under 
a Markovian approximation which assumes that the cou- 
pling of the cavity to both these systems is weak and that 
correlation times within these baths are much shorter 
than all relevant timescales, the density matrix of the 
cavity mode, p, evolves according to a quantum master 
equation (QME): 

P-W(p), (1) 
where the super-operator W [5] is of the Lindblad form 



[3IIS], 

W(p) = [l^pLI - \{LlL^,p}^ , (2) 
and the Lindblad operators read, 




= ^/v+la, (5) 
Li = ^/^ a\ (6) 



Here a, are the raising/lowering operators of the cavity 
mode, iVcx is the atom beam rate in units of the cavity 
lifetime (the inverse of the thermal relaxation rate) which 
we have set to one for simplicity, and v is the average 
thermal photon occupation number. The "accumulated 
Rabi angle" <^ (the effective vacuum Rabi frequency of 
the atom-cavity system multiplied by the time spent by 
the atom in the cavity [5]) quantifies the atom-cavity in- 
teraction. Li and L2 are the quantum jump operators on 
the cavity due to the interaction with the atoms describ- 
ing absorption of a quantum from an atom and passage 
of an atom through the cavity without absorption respec- 
tively. The other two jump operators, L3 and L4, are due 
to interaction with the thermal bath. 
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A. Steady-state properties 



The stationary state density matrix of the QME (l][6) 
is easy to obtain exactly. It is diagonal in the number 
basis, and reads [H |S] : 



Ps.s.{n) = Ps.s.(O) II 



m—1 



iVox sin^ i(t)y/m) 



1 



(7) 

An example of the steady-state behavior of the cav- 
ity as a function of the reduced atom pump parameter 
a — ipy/Nex is shown in Figs, [sjb) and (c). As the 
pump parameter is increased from zero, the average num- 
ber of cavity quanta, (n) , displays a series of local max- 
ima (minima) corresponding to situations where the cou- 
pling between atom and cavity is such that the atom ap- 
proximately performs a half (complete) Rabi cycle when 
traversing the cavity. As the effective Rabi frequency de- 
pends on the occupation of the cavity, the sharpness of 
this effect is smoothed out. 

The regions of sudden increase in (n) that occur 
around a ~ 1 and a ~ 27r,47r, ... correspond to dy- 
namical crossovers in the cavity state and are marked by 
sharp peaks in the corresponding Fano factor (defined as 
{{-n?) — {n)"^) / {n)) . The crossovers at a ~ 1 is classified as 
continuous and the others as first-order based on an anal- 
ysis of the behavior of — In [ps.s.('T')//Os.s. (0)] in the spirit 
of Landau theory [2| and on the observation that the re- 
gions of sudden increase in {n) around a ~ 27r, 47r, ... (but 
not the one at a ~ 1) get progressively sharper as N^x is 
made larger. Thus the crossovers are usually interpreted 
as being analogous to thermodynamic phase transitions 
in the limit N^x — > oo [31 13]. The number of minimas 
in — In [ps.s. [n) /ps.s. (0)] increases progressively as a gets 
larger, but each individual minimum typically becomes 
less distinct which leads to the progressive broadening of 
the dynamical crossovers [5]. 



B. Quantum Trajectories 

Instead of looking at the steady-state properties of the 
cavity density operator we will analyze the statistics of its 
quantum trajectories. The quantum jump trajectories of 
the system correspond to the time record of projection 
events due to the action of the Lindblad operators ^ 
|6]). We wish to classify such trajectories according to 
the number of these events |8]. The probability Pt(K) 
to observe K events after time t is given by Pt{K) = 
Tr[p'^)(i)], where p^^\t) is a reduced density matrix 
obtained by the projection of the full density matrix (i.e. 
cavity and baths) onto the subspace of K events [16] and 
pit) — X^if P^^"* (0- Here we will consider in particular 
quantum jumps associated with the jump operator Li 
Eq. ([3]) which describes the detection of unexcited atoms 
emerging from the cavity, as this type of measurement is 
readily performed in experiment. However, in principle 
we could define a counting variable associated with any of 



the four quantum jump operators and the results would 
be equivalent. 

For large times Pt{K) acquires a LD form |10) : 



Pt(i^)=Tr 



,(^)(t) 



-tip{K/t) 



(8) 



The function ipiKjt) (k = K/t) is a "large-deviation" 
function. It encodes all information about the probability 
of K at long times [10]. An alternative description for 
the statistics of K is provided by the generating function. 
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(9) 



The "counting field" s is the conjugate field to the observ- 
able K. In this dynamical context K and s are what pres- 
sure and volume, or magnetization and magnetic field, 
are in thermodynamical contexts. Although only the 
quantum trajectories where s — are observed in an 
experiment, we can nevertheless construct a generalized 
family of quantum trajectories, where the probability of 
observing K events after time t, P^^\k), is biased by 
the value of s [SI Hi], 



(10) 



In the long time limit the generating function also ac- 
quires a LD form [TU] , 



(11) 



The LD functions (p{k) and 9{s) are to trajectories 
[H [13 [T3] what entropy density and free-energy density 
are to configurations in equilibrium statistical mechanics 
[11) . The two LD functions are related by a Legendre 
transform, 9{s) — — min^ [(p{k) + ks] [TOIIIH]- The num- 
ber of events K is a time-extensive quantity. The scaled 
activity k = K/t is a dynamical order parameter as it 
serves to qualify and distinguish dynamical phases. 

Our approach differs from that often taken with stud- 
ies of counting statistics |18j as we focus on exploring the 
behavior beyond the s = limit. The function 0(s) has 
the convexity properties of (minus) a free-energy. Fur- 
thermore, its analytic properties as a function of s en- 
code non-trivial fluctuation properties of dynamical tra- 
jectories. In particular, singularities in 0{s) correspond 
to dynamical (or space-time [HI [H]) phase transitions. 
Just like in equilibrium statistical mechanics, our aim is 
to compute this free-energy-like function to uncover the 
dynamical phase structure of the micromaser. 

The statistical properties of the s-biased trajecto- 
ries are described by the master equation for Ps{t) = 



(12) 
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where the generahzed quantum master operator 
(GQMO) reads [HmS]: 

1 

+ Y.L,p,Ll--Y^{LlL,,p,}, (13) 



The LD function 0{s) is given by the largest eigenvalue 
of Ws 8J and it is then straightforward to obtain the 
moments of the distribution P^^\k). 



III. MEAN-FIELD ANALYSIS 

Although a full analytic solution for the steady-state 
density operator is readily obtained for s — this ap- 
proach is not easily generalized to the case s ^ 0. We 
therefore proceed by using a mean-field approach to cap- 
ture the behavior of the LD function. We check the accu- 
racy of this calculation using numerical methods in Sec- 
tion HVl 



The LD function can be estimated from ( 13 ) through a 



variational ansatz. For simplicity we start by considering 
the limit of zero temperature, = 0, and we restrict 
the analysis to density matrices which are diagonal in 
the number basis. With the rescalings a — > a^/W^ and 



the GQMO reduces to an operator, 



,t V )_ 



sin^ {a\/ a^a + — a) a 



+ \J flta + 8 a 
(14) 



where S = N^^. 

A variational approximation to the largest eigenvalue 
of Ws is obtained by maximizing Wg w.r.t. a and a\ 
that is, by maximizing the matrix element of Wg be- 
tween variational coherent states. We set a — e^^\/n 



and = 



/n, and the corresponding Euler-Lagrange 



equations are dWg/d^ = and dWg/dn — 0. The first 
of these can be solved, 



dWs 



s/2 fZ. |sm("%AW^)| 



= 



a) = e''l'^\Jn{n + 8) |csc {a.\/n + | 

(15) 

By replacing these relations into we obtain a varia- 
tional "free-energy" 



Tsin) = n — 2e "^^Vn sin (^aVn + S 



whose minimum w.r.t. n gives the variational estimate of 
the LD function. 



A similar procedure can be applied at finite temper- 
ature, V 0. In this case the variational free-energy 
reads. 



X + 6)v + e^^ siv? [a\fn + 5^ 

+ {n + 6)v + siv? {a^/n + 5^ . (18) 

The LD function obtained from the variational free en- 
ergy has multiple singularities. Fig. 3(a) shows the aver- 
age photon occupation of the cavity, n*, that minimizes 
the variational free energy (18), F'g{n^,) = 0. The left 
hand side of the plot, s < 0, displays multiple transitions 
where changes discontinuously. On the right hand 
side, s > 0, there is a single first-order line. In contrast 
to the s < transitions, this line ends at a critical point, 
Sc ^ and ac > 1. The existence of this critical point is 
most easily established from the form of at zero tem- 
perature, Eq. ( 16 ). In the limit of 5 = and at s = this 



variational free energy reads, J-s{n) {\/n — |sin {ay/n)\) . 
For a < 1 this function is minimized by — 0; for 
a > 1 it also has a second minimum at n* = sin^ (a^/n). 
At the critical point ^ \ these two extrema coalesce. 
Furthermore, in the limit of 5 = the coexistence line is 
along the s = axis. When 5 the critical point shifts 
to Sc > and etc 1 and the coexistence line bends to 
the right. 

The phase structure revealed by is that of the dy- 
namical phases. The LD function d{s) is the moment gen- 
erating function for and {k){s) = {K){s)/t = —0'{s). 
From Eqs. ( |17|18p we obtain the relation between and 

= (k) (s) /Ncy^ in the mean-field approximation: 

^ ^ e-'^^^n,{iy + 1) sin' (gy/^^TT?) ^^^^ 

(n* + 5)v + e^^ sin' [a^/n^ + 5) 

The activity /c* has the same phase behavior as as 
shown in Fig. 3(b). The close similarity between the two 
is easy to understand in terms of the energy balance in 
the system: at low temperatures the number of quanta 
absorbed per unit time by the cavity is dominated by 
the contribution from the atoms, (fc), whereas the rate of 
energy loss (to the bath) is proportional to (n) . 

For a > TT and s ~ the phase-diagrams of Fig. 3 show 
that the first-order transition lines converge to s = with 
increasing a. The mean-field solution of (18 1 predicts 
that these transition lines all accumulate at s = 0. 

The mean-field treatment of this section will become 
less accurate for larger values of a where non-linearities in 
the operator (14) become more prominent. The approx- 



(16) imation above assumes, firstly, that averages of products 
of a and can be reduced to products of their averages. 



and, secondly, that the operator ( 14 1 is normal ordered 



9[s) ~ — minJ^s(n). 



(17) 



Furthermore, the approximation assumes that the eigen- 
state associated to the maximal eigenvalue 9[s) is a co- 
herent state. This means that when a 3> 1 the mean-field 
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FIG. 3: (Color online) (a) Space-time phase diagram from the mean-field approximation (N^x ~ 100 and f — 0.15). The 
average cavity occupation (n) /A^ex is plotted as a function of a and s. (b) The same plot but for the (scaled) activity per unit 
time, {k}/Ncx- (c) Comparison of exact (black) and mean-field (red) phase boundaries. 



treatment may not give accurate estimates of (n) , for ex- 
ample near the dynamical crossovers at a ~ 2tt, Att, ... at 
s = [see Fig. 2(b)]. Nevertheless, the dynamical Lan- 
dau free-energy ( 18 ) does capture the singular changes in 



0{s), so we expect the mean-field treatment here to give a 
good approximation to the dynamical phase-behaviour of 
the micromaser, as is confirmed in the next section where 
we compare to the exact numerical diagonalisation of the 
GQMO. 



IV. EXACT NUMERICAL DIAGONALISATION 
OF THE GQMO 

For the numerical derivation of the LD function we 
constructed the matrix representation of the GQMO for 
given A'ex and v in the number state basis. The matrix 
was truncated at a state with occupation number N^^^x 
significantly larger than iVcx and the largest eigenvalue 
was obtained as a function of s and the pump parameter 
a. The activity K and other quantities were calculated 
by taking numerical derivatives of the LD function. 

The resulting phase diagram is shown in Fig. [T] There 
we show the average activity as a function of pump pa- 
rameter a and counting field s. The dynamical or space- 
time phase behavior is very similar to that obtained from 
the mean-field approximation. Fig.|3][c) shows a compar- 
ison of the exact and mean-field phase boundaries. These 
are very close to each other, and only start to disagree at 
large a: this is to be expected since the larger a the more 
prevalent non-linearities become which are captured less 
accurately by the mean-field expansion. 

Figure [3] shows the scaled average cavity occupation, 
(n)/iVox, as a function of pump parameter for different 
s. At s = the curve is smooth; this reproduces the plot 
shown in Fig.[2](b): the photon occupation oscillates with 
varying pumping, in a way which is suggestive of phase 
transitions without actually being singular at finite Alex- 
in contrast, for s ^ we find discontinuous behavior even 
for the relatively small value of N^x = 100 chosen. 

For s < the curves display singularities at the values 




FIG. 4: (Color online) Average cavity occupation, n (scaled 
by Ncx), as a function of pump parameter, a, for different 
values of the counting field s (see Figs. 1 and 3 for refer- 
ence). (Top) Active side of the dynamics (s < 0): n is dis- 
continuous at the first-order transitions (their location indi- 
cated by dashed lines). (Bottom) Inactive side (s > 0): here 
s > Sc ~ .0087, and there is a single first-order jump in n. 
(Center) For comparison, the same curve at normal dynamical 
conditions, s — 0, displays no singularities. 



of a where the transition lines are crossed, see Fig.|3] At 
these points the average occupation and the dynamical 
activity jump by a finite amount, so they are first-order 
transitions. These first-order phase boundaries do not 
seem to end, see Fig. 1, so (n) and {k) show qualita- 
tively the same behavior for all negative s. The sizes of 
the jumps (n) at the first-order transitions get smaller 
at larger a which is reflected in the behavior at s = 
where the crossovers get more and more smeared out as 
a is increased. For s positive and larger than the critical 
value Sc, the occupation vs. pumping curve shows one 
singularity when the rightmost transition line is crossed. 

The critical point is considered in more detail in Fig. 
5] Panel (a) shows the dynamical activity as a function 
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FIG. 5: (Color online) Diverging correlation time at the criti- 
cal point, (a) Average dynamical activity as a function of s for 
three values of the pump parameter a approaching the crit- 
ical value Oc ~ 1.342. (b) Correlation time r obtained from 
the difference between the leading and subleading eigenvalues 
of the GQMO. (c) The correlation time appears to diverge 
at {sc,ac) as r ~ |s — Sc|~^, approaching both from the left 
and right (the left branch has been shifted down by a con- 
stant factor for clarity). Results from numerically exact di- 
agonalisation of the GQMO, Eq. ( |13[ ); the Ifilbert space was 
truncated at Nmax ~ 350, where A'^ox = 100. 



of s for subcritical values of a approaching ac ~ 1.342. 
Panel (b) shows the corresponding correlation time r for 
quantum trajectories. This is calculated from the gap 
between the leading eigenvalue of the GQMO, i.e. the LD 
function 9, and the subleading one 92- t = {9 — 6*2)^^. 
This is just like calculating the correlation length in a 
thermodynamic problem where the partition sum can be 
obtained from a transfer matrix |11) . The correlation 
time appears to diverge at the critical point as r ~ |s — 
Sc|"\ see Fig.jHJc). 



V. CONCLUSIONS AND DISCUSSION 

The approach we have taken in this paper offers a new 
perspective on the dynamical transitions which occur in 
open quantum systems such as the micromaser. Instead 
of focusing on the steady state solution of the master 
equation, we have examined the statistical properties of 
the quantum trajectories corresponding to the measure- 
ment of the atoms after they leave the cavity. We used 
the large deviation limit as a way of exploring the analytic 
properties of the moment generating function at non-zero 
values of the counting field, s. Surprisingly, we were able 
to identify discontinuities in the large-deviation function 



for non-zero values of s which we interpret as correspond- 
ing to first order phase transitions. We also found a crit- 
ical point where a first order transition line terminates 
and which we found displays the standard scaling behav- 
ior associated with a continuous thermodynamic phase 
transition. 

The non-linear character of the quantum master equa- 
tion ([1]) gives rise both to the phase transitions we find 
at s 7^ and the well-known dynamical crossovers which 
are found from an analysis of the steady-state proper- 
ties of the micromaser. For a ~ 1, one can think of the 
s = behavior as being controlled by the nearby criti- 
cal point at (ac,Sc)- Furthermore, we showed that the 
critical point moves towards the s — line as Nex is 
increased, reaching it in the limit N^^^ ~^ 00, which is en- 
tirely consistent with the interpretation of the dynamical 
crossovers found from the steady-state density operator 
as true phase transitions in this limit. The proximity of 
the first-order transition lines to s = at higher values 
of a will also have observable consequences in terms of 
dynamical fluctuations (think for example of the promi- 
nence of vapor bubbles in a liquid at conditions near 
liquid- vapor coexistence) . 

Even though it is just the unbiased trajectories that 
can be measured in an experiment, our predictions about 
the presence of singular behavior at non-zero s could still 
be tested quite easily as the statistics for all the biased 
trajectories can be generated when those for the unbiased 
case are known. However, a detailed comparison with 
experimental data would require us to extend our model 
to include the imperfections that typically occur in real 
measurements of the atomic state. It may also require 
the issue of finite-time corrections to be considered. 

In this paper we have restricted ourselves to studying 
the simplest model of a micromaser and to examining 
only the most basic statistical properties of the trajec- 
tories. It will be interesting to examine more complex 
models where quantum coherences play an important role 
and the dynamics is not necessarily Markovian, as well 
as examining the dependence on s of higher cumulants 
of the counting distribution. We expect that the real 
utility of the method will be in analyzing more complex 
systems where a full understanding of the steady-state 
density matrix is not easily obtained. 
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